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I. INTRODUCTION 

The distance-redshift relation for far away objects 
plays an important role in cosmology. It has led Hub- 
ble, or rather Lemaitre to discover the expansion of 
the Universe; and the distance-redshift relation to far 
away Supernovae type la is at the origin of last year's 
Nobel Prize in physics for the discovery of the acceler- 
ated expansion of the Universe 0, Q . 

A next step that has been initiated recently considers 
the angular and redshift fluctuations of the luminosity 
distance, which may also contain important information 
about our Universe [3-0] • One important unsolved prob- 
lem is the question how strongly the distance-redshift 
relation may be affected by the fact that the actual Uni- 
verse is not homogeneous and isotropic, but the matter 
distribution and also the geometry have fluctuations. To 
first order in perturbation theory these fluctuations can 
average out in the mean and are therefore expected to be 
small. 

However, it has been found that they arc significantly 
larger than the naively expected value that would be of 
the order of the gravitational potential, namely, ~ 10~ 5 . 
An analysis in first order gave fluctuations of the order of 
10~ 3 , hence 100 times larger than the naive estimate 
Recently, Ben-Dayan et al. || have calculated a second 
order contribution to the distance-redshift relation of the 
order of ~ 10 -3 . Evidently, if the second order term is 
as large as the first order, this means that perturbation 
theory cannot be trusted. On the other hand, fully non- 
linear toy models, which have been studied in the past, 
always gave relatively small modifications of the lumi- 
nosity distance if the size of the fluctuations, spherical 
voids Q or parallel walls [ioj . is small compared to the 
Hubble scale. Hence the problem remains open. 

So far, the perturbative analyses of the distance- 
redshift relation have concentrated on scalar perturba- 
tions. In this work, we want to study the contributions 
from vector and tensor perturbations on a Friedmann- 
Lemaitre (FL) universe. This is interesting for several 
reasons. First of all, tensor perturbations are generically 
produced during inflation, and hence their contribution 
has to be added for completeness. Second, a passing grav- 
itational wave from some arbitrary source does generate 
a tensor perturbation in the distance-redshift relation 



to any far away object and could, at least in principle, 
be detected in this way. For single binary sources we 
have found that this effect is very small [Tl[; however, 
a stochastic background might lead to a detectable ef- 
fect. Even though vector perturbations are usually not 
generated during inflation (and if they are they decay 
during the subsequent radiation dominated phase), they 
are relevant in many models with sources like, e.g., cos- 
mic strings or primordial magnetic fields. A third impor- 
tant motivation to study vector and tensor contributions 
comes from the fact that at second order in perturbation 
theory, scalars also generate vector and tensor pertur- 
bations [H, HH . In a complete second order treatment 
these have to be included. With the formalism developed 
in this work, such an inclusion is straight forward. We 
plan to report on the result of these second order contri- 
butions in a forthcoming paper [l4| . A similar program 
is carried out in Refs. [IjjIlSl- There the authors discuss 
scalar, vector, and tensor perturbations and split them 
into E and B modes. The treatment of these papers is, 
however, more adapted to describe distortions of surveys 
and weak lensing, but the convergence calculated there 
is related to our distance fluctuations. 

The paper is organized as follows. In the next section 
we discuss the luminosity-rcdshift relation perturbatively 
at first order. In Sec. |TTT] we apply these results to ten- 
sor perturbations. We first derive the general first order 
expressions, which we then expand in spherical harmon- 
ics. We also give a numerical example for the gravita- 
tional wave background from inflation. In Sec. IIVI we 
treat vector perturbations and in Sec. [V] we conclude. 
Some lengthy calculations and some details are deferred 
to four Appendices. 

Notation: We use the metric signature (— , +, +, +). We 
denote the derivative w.r.t. the conformal time 77 with a 
dot. 



II. THE DISTANCE-REDSHIFT RELATION 

For an arbitrary geometry, defined through the metric 
g, a distance measure D from a source moving with 4- 
velocity us = u(xs) and an observer moving with 4- 
velocity uq — u(xq) can be obtained as a solution of the 



2 



Sachs focusing equation [l7j : 
d 2 D 



dX 2 



D. 



(1) 



Here A is the affine parameter of a lightlikc geodesic x fl (X) 
from the source to the observer, x^(Xs) = x^, x^(Xo) 
'a- 



x 1 ^, and 



1 



with fc M 



~dA~ ; 



(2) 



is the 4- velocity of the lightlike geodesic, and £ is the 
complex shear of the 'screen' defined below. The source 
and observer are made out of baryons; hence we identify 
the 4- velocity field u M (x) with the (baryonic and dark) 
matter velocity field. 

Considering a thin light bundle with vertex at the 
source, the luminosity distance is given by 



D L = (l + z)D, 
where z denotes the source redshift, defined by 



1 



k^u u \ Q loq 



(3) 



(4) 



We are considering past light cones without caustics 
between the observer and source positions. This is well 
justified as we are treating small perturbations on a 
Friedmann background. See Ref. [18[ for more details 
on the effect of caustics in the past light cone. 

The complex shear of the light ray bundle, S, is defined 
as follows [ijj : Consider two spatial orthonormal vectors 
e\ and which are normal to both the 4- velocity uo and 
k at the observer position and which are parallel trans- 
ported along k, such that Vk&a = for a = 1,2. The 
vectors e\, e-i are a basis of the so-called screen. Note 
that we do not require that u be parallel transported 
along k; hence e%, ei are in general not normal to u else- 
where than at the observer. The complex shear is defined 
by 



g (V e e, k) , with e = e\+ie2- 



£ = ig(e,V £ fc) = -i 

(5) 

We consider a light bundle with vertex at the source 1 . 
This leads to the following initial conditions (more de- 
tails are found in Appendix [5} for the Sachs focusing 
equation ((T|) 



D{X s ) = 0, D'(X s )=oJs = -9^k> l u' i \, 



(0) 



In a perturbed FL metric the Sachs focusing equa- 
tion ([!]) reduces, at first order, to 



d 2 D 
dX 2 



-1ZD. 



(7) 



1 A light bundle with vertex at the observer yields the angular 
diameter distance, which is related by a factor (1 + z) 2 to the 
luminosity distance which we determine here. 



Since the complex scalar shear £ vanishes for a confor- 
mally flat spacctime, |S| 2 contributes only at second or- 
der. To first order in 1Z, Eq. Q with initial conditions (|6]) 
is solved by 



D{Xp) 



= (A -A s )- { °dxf dX"Jl(X-X s ) 

J As J As 

= (Ao-A s )- / °dX(X-Xs)(X -X)K, (8) 



where we have used the identity 

Vo rV rno 

dr, / d V 'f (r/) = / dr, ( Vo - rj) f (rj) (9) 

r/s J VS J VS 

for the second equal sign. 

Of course, in a perturbed FL universe 1Z is not first 
order; it also has a zeroth order contribution. But a 
perturbed FL universe is conformally related by the scale 
factor a to a perturbed Minkowski spacetime and lightlike 
geodesies are invariant under conformal transformations. 
Two conformally related metrics, 



have the same lightlike geodesic curves, and only the 
affine parameter changes, dX — a 2 dX, such that k* 1 = 
a~ 2 k^. Also the (normalized) matter 4-velocity changes, 



a u> 



so that the redshifts are related by 



z + 1 = — (Sz + 1) , where — 
as a s 



z+l 



(10) 



is the background redshift, i.e., the redshift in an un- 
perturbed Friedmann-Lemaitre universe, and 5z is the 
source redshift according to definition ((4J w.r.t. to the 
perturbed Minkowski metric g, while S is the one w.r.t. 
to the perturbed FL metric g. We remark that z is the 
true (observed) redshift. In what follows we shall nor- 
malize the scale factor at the observer to one, ao = 1- 
The distance D is not affected by a conformal factor, so 
that the effect of the expansion on the distance simply 
leads to a rescaling Q 



D L = (l + z)D L . 



(11) 



We now compute the luminosity distance in a per- 
turbed Minkowski spacetime, Dl, and then relate it to 
the one in a FL spacetime, Z?£, by the above rescaling. 
Let (l,n l ) be the 0-order term of the lightlike velocity 
vector fc M (in the noncxpanding Minkowski spacetime). 
The lightlike condition implies |n| =1. We normalize 
the affine parameter A such that ljs = kg = 1 at 0th 
order. To determine the redshift Sz, we have to solve 
the perturbed geodesic equation for \i = only (in order 
to determine k to first order), since the peculiar veloc- 
ities are already first order. The Christoffel symbols of 
Minkowski space vanish, so that the geodesic equation 
for fi = to first order is simply 



dk° 
~dX 



1 oo 



2r' 



,0' 



r?,nW 



0. 



(12) 
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We normalize the affine parameter A such that k s = 1, 
and Eq. ([T2")) is solved by 



k° 

Kq 



dX(T° 00 + 2T i0 n i + T ij n i n j ). (13) 



The geodesic equation (|T2j) will be useful also in order to 
express the distance D in terms of the conformal time r\ 
instead of the affine parameter A. For this we use 



1 



dA'(r° + 2rV l + r^V), (u) 



A,. 



which, in first order, leads to 



Xo-Xs = Vo^lS+ d V d V '{Tlu + 2Y%rt+Y%n l n3) . (15) 



The conformal time and the background redshift are 
not observable. We want to write the distance as a func- 
tion of the true (observed) redshift z = z + 5z, where 
Sz = (1 + z) Sz according to Eq. ([TU|). Following the ap- 
proach presented in Q we compute 



D L (rj s ,n) = D L (r? (z) , n) = D L (z, n) 
= D L (z,n) - -^zD L (z, n) 



JS, (16) 



with 

-f=D L (~z,n) 
dz 



where Hs = — 



= — Dl (z, n) + first order 

az 

= (1 + -D^ + "Hg 1 + first order, 

(17) 



In other words, we evaluate the distance at the true (ob- 
served) redshift Dl (z, n) by using Eqs. (fT6llT7|) in order 
to relate £>i(z,n) to £>L(?7s,n). 

From Sec. IIIII on, to simplify the notation, we denote 
the true (observed) redshift with z instead of z. We shall 
not use z anymore. 



III. THE DISTANCE— REDSHIFT RELATION 
FROM TENSOR PERTURBATIONS 

We first consider a perturbed Minkowski metric with 
tensor perturbations only, defined by 



ds 2 



-drf + (Sij + 2Hij) dx l dx\ 



(18) 



where the tensor perturbations arc divergence-free 
Hj i = 0, traceless H\ = 0, symmetric Hij = Hji, 
and spatial H^q = 0. By definition, a spin-2 pertur- 
bation is gauge-invariant. To use a notation consistent 
with the next section, we introduce the gauge invariant 
shear on the {t = constant} hypersurfaces er^ = Hij 
(see, e.g., HJ). 



A. The perturbation equations 

From the Ricci tensor calculated in Appendix [B] wc 
obtain 

Tl = --n i n j nH ij , where \3 = d' t d li . (19) 

Note that this is the Minkowski space d'Alembertian, 
without expansion. The geodesic equation (|13[) for /i = 
leads to (see Appendix [B] for details) 



k = 1 — / dX aijU 1 ^ . 



(20) 



We consider the 4- velocity (u^) = (1, 0) because the spin- 
2 perturbations can not source peculiar velocities at lin- 
ear order, so that we obtain the redshift to first order 

1 f Xo ■ 

1 + Sz = w = l+ dXaijrtn 3 . (21) 

With Eqs. ((SJ [TTJ [T5|) we find the luminosity distance in 
a perturbed FL universe with ds 2 = a 2 ds 2 , as a function 
of the background redshift 

Dl (z, n) = (1 + z) (i] - f]s) 

x ( 1 + / drj (7ijn % n^ 
'ns 

r'/o 



, rjq - T} 
drj <JijU n 



>,,l 



vs Vo ~ VS 

v ° dy^-^^-^n 

vs Vo ~ Vs 



(22) 



We have again used (O to reduce the double integral. 
Wc finally express the luminosity distance in terms of the 
true, observed redshift z. Using Eqs. (|16I17[) . wc obtain 

D L (z,n) = (l + z){r)o-Vs) 



x 1 



'/<> 



Vo -Vs j vs 

f)o 



drj (Jijtfn 3 



, vo-v i j 

drj <7j 7 -n n J 

s Vo-Vs 



'is 



Vo - Vs 



(23) 



The origin of the different terms in the redshift-distance 
relation is as follows: the first line is the unperturbed 
expression for the luminosity distance in a FL universe 
at the observed redshift z, the term on the second line 
derives from the redshift correction, the one on the third 
line from the relation between the conformal time 77 and 
the affine parameter A, and the one on the last line from 
the Sachs focusing equation. We can interpret this last 
term as a lensing effect. The first two terms come from 
the perturbation of the redshift. 
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For a fluid with a vanishing anisotropic stress the 
rcdshift-distancc relation becomes 

D L (z, n) = (1 + z) (r]o - Vs) 

x ( 1 2 / dr] Oijn n 3 



Vo - Vs J ns 
no 



+ / drj 



vo-v i 3 



n s Vo - Vs 
n ° d ^-Vs)(vo-v) n ninjai . 



> vs Vo - Vs 

where we used the Einstein equation [2 



H i:i + 2UH i0 - \7 z H l3 = 0, 



and Eq. (|T9|) to replace TZ. If the cosmic fluid is not ideal, 
but has anisotropic stresses, these add to the right-hand 
side of Eq. ([2"5"]) (see [20| ) and correspondingly to the last 
line in Eq. ([2~4")) . the lcnsing term. 



B. Spherical harmonic analysis 

We want to determine the power spectrum of the lumi- 
nosity distance. In the unperturbed FL background the 
luminosity distance to the redshift z is given by 



D L (z) = (l + z) (vo - Vs) ■ 



(26) 



We define the relative difference in the luminosity dis- 
tance as 



A L (z, n) 



D L (z,n)-D L (z) 



Dl (z) 

1 f7° 



u 



VO - VS Jr,s 

+ (r]-Vs)(vo-r])'H 



s 1 + (vo-v) + 
a.jn'n' . (27) 



Note that we evaluate the unperturbed distance at the 
true, observable redshift. 

We are interested in the angular power spectrum of this 
observable, q (z, z'), which depends on the redshift of the 
two sources and is defined by the two point correlation 
function 

(A L (z,n)A L (z',n')> 

= ^E(^ +1 ) c ^ z ' z, ) p ^ n - n ')- ( 28 ) 

£ 

In the distance-redshift relation ([231) [and, in particu- 
lar, in Eq. ([2~4"]) for an ideal fluid] we have several times 
the term tfnPtJij (77, x (77)) where x (rf) = xo-n (770 — v)- 
In terms of its Fourier transform this is 



n % n>Oij (t7,x(t7)) = 



d 3 k 



o-yfokJnVe - * 1 ™^. (29) 



Without loss of generality we choose xo = 0. Setting 

k = k |k| = kfc, fi = k ■ n, A77 = 770 - V, (30) 
we obtain 

d 3 k 



n l n J aij (rj,x(r])) 



(27T) 



j^^kJnVe^. (31) 



Writing the exponential in terms of spherical Bessel func- 
tions 



D ifj.kAr) 



£(2* + l)i< u (k At,) Pt(fi), (32) 



1=0 



(25) we find 



d 3 k 



Way + 1) i e h (kA V ) P e (fi) . 

(33) 



With respect to a helicity basis in Fourier space 

{e(+),e(-),k} , eW = i=(e 1 ±ze 2 ) , 
such that 



e(±) • n : 



1-/J 2 



we have 



an e* <g> e 3 = <t+ e (+) ® e (+) + a e (-) <g> e ( ~ 



(34) 



(35) 



(36) 



We introduce the spherical harmonics with respect to 
some arbitrary z direction given by a unit vector e as 
Y£ m (n, e), since we shall use them w.r.t. different z axes. 
The addition theorem of spherical harmonics is 



Pt(») 



-in 



2i- 



r ^F/ m (k,e)y £m (n,e). (37) 



Using the following spherical harmonics definition 



2tt 2 



, (38) 



Y2 ±2 ( n,k ) = \ / — sin 
' V 'J V87T 

we can rewrite Eq. ([33)) as 

n l n J a i: j = J -^-yg (<5- + F 2 2 (n, k) + <x~F 2 - 2 (n, k 

x V % E 4?r/ ^ ( fcA7 ?) Y tm (k, e) y £m (n, e) . (39) 

We now introduce the initial tensor power spectrum 
Ph (k) through 

(H± (ns^H^^y)) 

= (2tt) 3 5& ( k - k') P H (k) T k (775) T k (jis>) , (40) 
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where T k (77) is the transfer function with the initial con- 
dition T k (-q) — >(fcr7— )-o) 1- The j( 3 )(k—k') is a consequence 
of stochastic homogeneity. For the shear we then have 

= (2^) 3 ( 5( 3 '(k-k')PH(fc)T fc (?7s)r fe (^0- (41) 

Next, we express the terms in the distance-redshift rela- 
tion with the help of the power spectrum of the integrand, 

(n i n j a ij n' l n' k ai k ) = ±- £ (21 + 1) 0,(77, rf)P t (n • n') . 



A lengthy but straight forward calculation yields |2C 



1 (£+2)\ 



,Jt(kArj) u(kArf) 



(kArj) 2 (kArj') 2 ' 
, (43) 

Using the Limber approximation (see Appendix [TJJ 
for the time integrals in Eq. (f2~4"|). and, in particular, 
Eqs. (jDSp . (|D6|) and (|D7|) . we can simplify the time inte- 
grals, and we find the coefficients (under the ideal fluid 
assumption) 

1 (1 + 2)! J? 1 1_ 

CHZ ' Zj " 7T (£ - 2)! £ 4 Arys Afj S ' 

/•GO 

x / dkP H (k)T 2 ( m , k ) (A + BH( Ve , k ) + C-H 2 ( m , k )), (44) 
Jk* 

where we have introduced 



A?75 =rio-VS 
Ve,k = Vo 
k* = 



I 

" k' 


2 1.58 


max< 


' £ I \ 


k Aj] S ' A?7s/ J ' 



At] S ' = Vo - T)g>, (45) 
(46) 

(47) 



and 



A = H s H s , + ^ - - (H s + 

B = 



fc 2 fc 



(48) 



+- (A ?7S + A ?7S , + n s 1 + Hg, 1 ) , (49) 

I 2 £ 3 £ 4 

C = —A Vs Ar ls ,- — (Ar 1s + Ar ls ,) + —. (50) 

More details can be found in Appendix |TJ] 



C. Application 

As an example, we consider a flat primordial tensor 
power spectrum as expected from inflation Pjj (fc) = 
a/k 3 . If r denotes the tensor to scalar ratio, the scalar 
amplitude as measured by the Wilkinson Microwave 



Anisotropy Probe experiment [22[ yields a ~ r x 10 9 , 
such that the tensor power spectrum becomes 



P H (fc) ~ - 10- 9 . 



(51) 



Considering an ideal fluid, the transfer function T k (77) is 
the solution of the differential equation (|2"5"|) . 



f k (n) + 2UT k ( V ) + k 2 T k (77) = 0, 



(52) 



with initial condition T k (rii n ) = 1 and T k (r]i n ) = for 
kr]i n <C 1. In a matter (or radiation) dominated universe 
this differential equation can be solved analytically in 
terms of Bessel functions. The growing (not decaying) 
mode is given by 

T k (v) = {kr,) 1/2 ~ q Y 1/2 _ q {kr 1 ), where a oc r? 9 , (53) 

and Y v is the Bessel function of the second kind of order v. 
At late times, when the cosmological constant dominates, 
we cannot write the scale factor a (rj) as a power law and 
we have no analytic solution to (|52")) . To determine the 
Q coefficients, we have solved the differential Eq. (|52")) 
numerically. 

The resulting power spectrum ci(z,z) for different 
source rcdshifts is shown in Fig. [TJ 
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FIG. 1: We show the tensor power spectrum rescaled by £ 4 for 
the fluctuations in the luminosity distance for different values 
of the source redshift (z = 0.5, dotted pink line; z = 1, dot- 
dashed blue line; z — 2, dashed green line; z = 3, long-dashed 
orange line; 2 = 4, solid red line). In the figure we have set 
r = 1. 

Clearly, for sufficiently large £, ct(z,z) oc £~ A . The 
simplest way to understand this scaling is to note that 
once a mode enters the horizon, the tensor fluctuations 
scale like J adrj ~ H cx a k /a oc (kr/)~ q , where a k = 
a(i] = 1/fc) denotes the value of the scale factor at horizon 
entry. For modes that enter during the radiation era 
q = 1, while for modes that enter during the matter 
era q = 2. Hence J adrj oc H oc H- m /k q is acquiring 
a factor k~ q with respect to the scale invariant initial 
spectrum. This leads to a red spectrum, fc 3 (J a) oc 
k~ 2q and £ 2 ci(z,z) oc £~ 2q . This spectrum turns from 



6 



ci oc £~ A for scales that enter the horizon in the radiation 
era to q oc £~ 6 for scales that enter the horizon in the 
matter era. For z = 4 this happens roughly at £ ~ 20. 
Of course, the transition is quite gradual. 

Comparing Fig. [1] with the results from scalar pertur- 
bations we see first that the tensor contribution is 
much smaller, nearly 8 orders of magnitude. We ob- 
tain £ 4 c e (z) - 5 x 10" 13 for z = 4 and t ~ 40 while 
scalar perturbations yield £ 2 ci(z) ~ 10 -5 for z = 4 and 

I ~ 100. Furthermore, despite also being proportional to 
the lensing term, it scales differently with £. This comes 
from the fact that the scalar lensing term is determined 
by the spectrum of k 2l $>, where VP is the scale invariant 
Bardcen potential, while for scales that enter during ra- 
diation dominated expansion, aijn 1 ^ is suppressed by a 
factor of 1/k. 

Interestingly the tensor signal is not monotonic in red- 
shift up to z ~ 2. It has a sharp minimum at z ~ 1.65. 
To illustrate this, we also plot ce(z,z) as a function of 
the source redshift for different values of £ in Fig. [2] 




z 



FIG. 2: We show the tensor power spectrum rescaled by £ 4 
for the fluctuations in the luminosity distance as a function of 
the source redshift for different values £ (I = 60, dotted blue 
line; £ = 40, dot-dashed green line; I = 20, dashed orange 
line; £ = 10, solid red line). Also here r — 1. 

The signal drops to at z c = 1.65. This comes from 
the fact that it is dominated by two terms with opposite 
sign. To see this, we also show the contributions from 
the three terms in the square bracket of (|2T|) individually 
in Fig. H 

If the source redshift is small, rjs ~ r/o, the first term 
oc —Tig 1 ~ —rjs/2 dominates, while for large redshifts, 
Vs *C Tjo, the second term oc (?yo — f]) dominates. If 
aijtfn 3 has a definite sign, the result inherits this sign for 
small redshifts and the opposite sign for large redshifts. 
The sign changes happens around rjs = rjo/2 correspond- 
ing to a redshift z c ~ 3. This is not expected to be very 
precise; in particular, we have neglected the time depen- 
dence of the transfer function in this argument. The 
more precise numerical evaluation gives z c ~ 1.65. In- 
terestingly this redshift is close to the maximum of the 
angular diameter distance Da(z) = (rjo — 7js)/(l + z). 




FIG. 3: We show the contributions from the different terms 
of Eq. (|27p to the tensor power spectrum rescaled by £ 4 for 
the fluctuations in the luminosity distance as a function of the 
source redshift for I = 40. (First term: dotted red line; second 
term: dot-dashed brown line; third term: dashed blue line; 
correlator of the first and second term: solid green line) . The 
third term is always subdominant. We plot the correlator 
between the first and the second term (solid line) with the 
opposite sign since it is negative. The second and the first 
terms have opposite signs, and they cross at z c ~ 1.65. 



The results shown in Figs. Q] to [3] have been calculated 
with the following cosmological parameters: h = 0.7, 
H X = 2997.9/i~ 1 Mpc, fl m fi 2 = 0.13, fl r h 2 = 4.17 x 
10~ 5 , and Q,a = 1 - fi m — ^r- 



IV. THE DISTANCE— REDSHIFT RELATION 
FROM VECTOR PERTURBATIONS 

We now consider vector perturbations. As for tensor 
perturbations, we can divide out the cosmic expansion 
for lightlike geodesies. Hence we can consider Minkowski 
space with purely vector perturbations. The metric is 
then given by 

ds 2 = -dr] 2 -2B i dx i dT]+{6i j + H id + H hi ) dx l dx 3 , (54) 

where the perturbations are divergence-free, 
= H*i= 0. Using Bij = B (iJ) and R l3 = H m , 
where ( ) denotes symmetrization, the shear on the con- 
stant time hypersurface is given by (20| ct^ = + ijy 
or in 3-vcctor notation <Xj = Bi + Hi, and cr^ = crnj). 
This quantity is gauge invariant [2(il |. 

A. The perturbation equations 

With the Ricci tensor calculated in Appendix IC II wc 
obtain for vector perturbations 

K = i (V 2 (<T in l ) + ciynV) . (55) 
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To determine the rcdshift we first evaluate the geodesic 
solution (fT3"]) with the Christoffcl symbols derived in Ap- 
pendix IC 11 

/•Ao 

k% = 1 - / dX g, ,„■„ '■. (56) 

Vector perturbations can have a nonvanishing peculiar 
velocity term. We define the observer 4- velocity (u^) = 
(l,B l + v l ). The peculiar velocity v l defined in this way 
is gauge invariant. It is the vorticity of the matter flow. 
We now obtain 



the procedure used so far is completely geometrical. If 
one is interested in general relativity (GR), then the two 
gauge-invariant variables v % and cr^ are not independent 
but related via Einstein's equations [2(J , 

V 2 cr, = -WirGa 2 Vi(p + p), (61) 

where p and p are the background density and pressure, 
respectively. 

B. Spherical harmonic analysis 



1+Sz = ,o +w = 1+ni t vh - v s)+ [ ° dX • 
-K + riiV j Xs 

(57) 

After a short calculation, using the results of Sec. [Til wc 
find the luminosity distance 

D L (z,n) = (1 + z) (r)o - Vs) 

( r° ■ ■ 

x 1 + niV l — 2mVg + drj Oijtfn 3 



1 



'/<> 



dr)(r]o - lijOiirtn 3 



Vo - Vs J vs 

n ^-Vs)(vo-v) 7Z 



Vo - Vs 



Since we are interested in expressing the luminosity dis- 
tance as a function of the true redshift, we have to evalu- 
ate Eq. (|58[) at z and subtract the correction term defined 
through Eqs. CGI and ([iT]). 



D L (z,n) = (1 + z) (vo ~ Vs) 

— — £ n t Vo - TliVg I 



Vo - Vs 



^s 1 
Vo ~ Vs 



H 



-i 



no 



' dr\ aijrf n 3 



Vo - Vs J vs 

1 [ r >° ■ 

/ dr](r)o - V^ij^nP 

Vo - Vs J V3 

n ° drl (n ~ Vs) (VQ - V) n 



'/.-■ 



Vo - Vs 



(59) 



A £ (*,n) r- 



Vo - Vs 



UiV - TliV s 1 



^s 1 

Vo - Vs 



1 drj aijn n J 



VO - VS Jrjs 
1 flo 



dv(vo - VjOijrfn 1 



VO - VS Jrjs 

V ~ Vs) (rjo - V) ^ 



VS 



Vo - Vs 



(60) 



This expression depends only on the gauge-invariant 
quantities v l and it should. We note also that 

we did not assume any gravitational theory yet. Indeed 



As for the tensor perturbations, we are interested in the 
term n l n?Oij. The main difference is that in the vector 



case we have <7y = ou^\ in real space and &i 
in Fourier space. This leads to 



ik 



i l n j aue^ kAri = -kn j a 



aijt -- n,v " r d(k&v) 

With this we can write Eq. ([551) as 



(62) 



d 3 k 

(2nf 



k ricTjY^ (M + 1) * l j'e(k^v) Pt(p) ■ 



1=0 



(63) 

With the helicity basis defined in Sec. IIII1 the addition 
theorem of the spherical harmonics ([3"T]) , and 



Yi.±i(n,k 

we obtain 



-,±i4 



_3_ 

^4tt 



1 -p 2 



(64) 



X47T1 



d 3 k 

(2tt 
4tt 
3" 



3 fc(y 1 , 1 (n, k) a + -yi > _i(n,k)ff 
£Y j' e (kr)Y* n (k,e)Y (m (n,e). (65) 



If we assume that vector perturbations have been gen- 
erated at some time in the past, we can define the vector 
power spectrum as for the tensor case as 

(<7± (vs,-k)^* (vs*,*!)) 

= (2tt) 3 5& (k - k') P a (k) T k ( Vs ) T k {ns>) -(66) 

If we do not want to consider the case of early generation, 
we simply have to replace P<j(k)Tk(vs)Tk(vs') by a time- 
dependent power spectrum, P a (k,rjs,vs')- The model 
under consideration (e.g., cosmic strings) then has to be 
used to determine this time-dependent power spectrum. 
If, however, vector perturbations evolve freely, we can 
then compute the shear power spectrum as for tensors, 

(nVo-ynV*^) = J" £ (2* + 1) HV, rjf)Pt (n ■ n') , 



(67) 
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finding (see Appendix IC 2 



be defined as 



ft: 



21 (I + 1) 

7T {21 + if 



dk fc 4 P CT (fc)T fc (r,)T k (rj') 



x [jU (fcr)il-i (kr')+j' e+1 (kr)j> +1 (kr')] .(68) 

As mentioned above, in more realistic scenarios, where 
vector perturbations are generated, e.g., via anisotropic 
stresses from topological defects or by second order per- 
turbations, we obtain a shear power spectrum of the form 
P a (fc, 77, 77'), which cannot be factorized into a random 
initial spectrum and a deterministic transfer function. 

In Appendix IC 2l we nevertheless, for sake of complete- 
ness, continue with expression (|68[) to derive the vector 
angular power spectrum for the luminosity distance fluc- 
tuations. We do not repeat the lengthy, complicated, and 
not very illuminating formulas here. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we have calculated the angular power 
spectrum of the linear vector and tensor fluctuations in 
the distance-redshift relation. For vector perturbations 
we have simply derived the formulas and for tensor fluc- 
tuations we have applied them to an initial spectrum of 
fluctuations from inflation. It is interesting to see that 
the tensor-distance fluctuation spectrum is not simply 
suppressed by a factor r as one might naively expect, but 
by about 8 orders of magnitude more. The reason for this 
is mainly that tensor fluctuations decay once they enter 
the horizon, while, on the contrary, scalar perturbations 
start growing. We therefore expect that the tensor sig- 
nal generated from scalar perturbations at second order 
dominates over the small first order signal. The calcu- 
lations of these second order contributions are left to a 
future project flij . 

We have also found that the tensor signal is not mono- 
tonically increasing with redshift as we would expect it 
from a pure lensing signal. This is due to the fact that the 
total signal is the sum of a redshift part, cx— Sz/(z + 1), 
and a lensing part. At redshift z c ~ 1.65, which is 
close to the redshift where the angular diameter distance 
Da = (vo — ?7s)/(l + z ) has a maximum, these terms can- 
cel, and at higher rcdshifts the rcdshift-tcrm dominates. 
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Appendix A: Sachs focusing equation 

Much of this work is based on the Sachs focusing equa- 
tion [T3|. It has been shown [23| that the distance can 



D = v/|detP| 



(Al) 



where V is the Jacobi matrix that satisfies the differential 
equation 



d 2 V 
d\ 2 



-K - Re (F) -Im (F) 
Im (F) -11 + Re (F) 



V. 



(A2) 



with F = ^R afJl p v e a l^ k^k v ' . The determinant of the Ja- 
cobi matrix T> describes the area of the thin light beam 
and its square root is therefore a distance, (|A1|) . if the 
affine parameter A is normalized such that ujs = 1- In 
general one can use the Sachs focusing equation also with 
a different affine parameter normalization. Indeed from 
the distance definition 



D 



dAc 



and the solid angle aberration [l7[ 

dn 



k^ 



we find, setting ui = 1, 



D 



'dAc 



dAc 



uiD. 



(A3) 



(A4) 



(A5) 



Since we are considering a light beam with a vertex at 
the source position, the initial conditions of the Sachs 
focusing equation are 



D(\ s ) = 0, 



dD (A) 



dX 



= 1, 



(A6) 



A=Ac 



if we normalize A such that ws = 1. The general initial 
conditions for an arbitrary affine parameter A, are given 
by ([6]) . Choosing wj = 1 + z one obtains the luminosity 
distance while lo§ = (1 + z) _1 gives the angular diameter 
distance. 



Appendix B: Details for tensor perturbations 

Here we write down the nonvanishing Christoffel sym- 
bols and the Ricci tensor for the metric (fT5|) . 



r° 

'j 



V (Huj + Ftij.i — Hij,i) = Hij, 



(Bl) 



r}o = 7 f h (Hko,j + Hjk,o — Hjo,k) — Hji, (B2) 



=> r^ = 0, 



= H 

rL = h. 



Huj — H 



Hij.i ? 



(B3) 
(B4) 



T l ■ = 0. 

%3 
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Ro 



Rio — 0, 



Rij — ij 



II,, u — —OH i: 



These components lead to 



n=\(R m 



2R M n l + Run 1 !!, 3 



2 3 



(B5) 
(B6) 



(B7) 



Appendix C: Details for vector perturbations 

1. Christoffel symbols and Ricci tensor 

Here we write down the nonvanishing Christoffel sym- 
bols and the Ricci tensor for the metric (El 



Too — Bi 



T}o = ~(dr,(Hi d +H jti )-B i!j+ B jti ), (C2) 



(CI) 

(C2) 
(C3) 



I" ^{lh, ■ ■ II,, ■ II,,) a,,. (C4) 

Y% = 11,., ■ II,, ]Tr° =o, (C5) 

i 

r %] = II,.,, T\ 3 = 0, (C6) 

r jk = H i,jk, (C7) 

r ifc = H i,ik =>■ T\ k = 0, (C8) 

(C9) 
(CIO) 



-Roo 


= o, 




1 




2 




1 


Rij 


2 



With these Ricci tensor components we can easily com- 
pute 



n = - {Roo + 2R io n i + R 4j n i n 3 ) 



(C12) 



2. The eg coefficients 

We first derive in detail Eq. We use the relation 



F 1)±1 (n,k) = A /^ Yl Y lm ,(n,e) Tl Y 1 * m ,(k,e). (C13) 

m'— — 1 



Here Tl Y 1 * m ,(k;, e) is the vector spherical harmonic. The 
general addition theorem for spin weighted spherical har- 
monics used in Eq. (|C13|) above can be found, e.g., in poj . 



Applying this to Eq. (|55|) . we obtain 
3 - 

im m — — 1 



n l n 3 a. 



'j 



J2 Y lm ,(n)Y em (n) 



x 



d A k 



kj' t (kAr,)Y; m (t) 



(27T) 



(C14) 



Here we omit the arbitrary unit vector e in the notation, 
and we have introduced the helicity basis, 



defined in Eq. (|34j) . In the special case with n = e we 
obtain 



W^ 1 

rl 3 k 

kf e (kA v )Y; (k) 

(2tt) 



where we have used 

Yi m (e, e) 



2^+1 

47T 



(C15) 



(C16) 



and Eq. (|C13p for n = e which yields, 



Ki,±i(e,k) = 



Y Yi m <(e,e) Ti r l * m ,(k,e 

rn' — — 1 

:i*io&e). 



(C17) 



Since the two point correlation function (|67|) depends on 
the angle n • n' only we can set n' = e without loss of 
generality. With this we find 



i l n 3 ' Oij^e 3 ' <Jij) 



A All 



x J dk fc^fcA^fcA?/)^ (k)T k ( v )T k ( v ') 

x J / m ,(k)y/ m (k)y Zo (k), (cis) 

where we have introduced 

/w(fc)= _ x ^(k) _ilio(k)+ ^(k) x Y 10 (k). (C19) 

Since the spherical harmonics form an orthogonal basis 
on S2 , we can expand the product of two of them again in 
terms of spherical harmonics using the Clebsch-Gordan 
coefficients [24| . In the case of Eq. (|C18|) we use 



£+i' 



(2£ + l)(2l+l) 

L=\P-<"\ ' v 



x (e,0,I,0\L,0)(£,m,lO\L,m)Y£ m (k). (C20) 
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The dependence of the spherical harmonics on the az- 
imuthal angle <fi, 



v* -y p -tm4 

1 Lm e 



and 



±1 Y* m , cx (C21) 



implies that the integral over angles in Eq. (|C18|) only 
contributes for to' = —to. Therefore also ms{— 1,0,1}. 
Since f m i (k) contains only terms that either do not de- 
pend on 6 or that are quadratic in sin (9) and cos(0) the 
only nonvanishing contributions arc L = or L = 2. 
Analogous to Eq. (|C20|) we write 



Yl %m (n) Yi _ m (n) = ^ 



3(2^+1) 



4?r(2n+ 1) 

x (£, 0, 1, 0|n, 0) (£, to, 1, -m|n, 0)Y ra0 (n) . (C22) 
The addition theorem for the spherical harmonics implies 



P n (n • e) 



4tt 



2n + 1 



(n) . 



(C23) 



Using these identities we can rewrite the correlation func- 
tion (|C18j) as 



n l n 3 tjije l e l ai 



X (21+ 1) [21+ 1 (2L + 1) 



EEEE 

n L— 0,2 ^ / m— — 1 
-1/2 



3tt 3 / 2 



x(€,0,I,0|L,0)(£,to,€,0|L,to) 
x (£, 0, 1, 0|n, 0) (£, m, 1, -m|n, 0) 

, dA : A : 4 j;(fcA ?? )4(fcA77')P CT (A : )r fe ( ?? )T fc (? 7 ')SLrnP n (n-e), 

(C24) 



where we have introduced 



B 



Lrn 



dn k /_ m (k)F£ m (k). 



(C25) 



The nonvanishing coefficients are given by 
5oo = ^, B2,±i = U^-, B 2fi = --^=. (C26) 

V7T 2 V 07T V57T 

Computing the sum of the Clebsch-Gordan coefficients, 
we find 

(^eV CT ,)=^^iip,(n-e) 
x / dfc k A P a {k)T k { n )T k { n ') 

x [^(fcA^jt^fcAr/) + J ; +1 (fcAr ? ) J ; +1 (fcAr/)].(C27) 

Using this result for the q's defined in Eq. (|57)) . we obtain 
directly Eq. (pgj) . 



Now wc compute the full a coefficient defined in (|2"5)l . We use the Limber approximation (Appendix [D]) to do the 
time integrals J dr]Ti i: ('r])j' e _ 1 (k(rio — v))- Wc start with the Doppler terms, the first line of Eq. (|60[) . which contribute 
to the dipole term only 

d _ 47r 'H s 1 H s } 2 / P9S n 

c i - ~q~ "T — -T < \vo\ > ■ 

3 A?7s A775/ 

This dipole is the same as the one from the scalar analysis We cannot, of course, decide which part of the observer 
velocity comes from scalar perturbations and which part from vector perturbations. Since this dipole term is highly 
nonlinear, we neglect it in the subsequent analysis and consider only £ > 2. We now determine the other terms. From 
the peculiar velocity of the source we obtain 



ji) _ 2£ {£+!) 



n(2£ + iy 



1 - 



ATys 



1 - 



A77S 



-) J dk k 2 P v (k) n ( vs ) T% (^0 Y, h ( fcA? ^) h ( fcA ^') . (C29) 



where we have introduced the velocity power spectrum defined by 

(« ± (775, k) o±* ( Vs > , k')> = (2tt) 3 5^ (k - k') p v (k) n ( Vs ) n (vs-) 

and ■0 ± is the peculiar velocity in terms of the helicity basis defined in Eq. (f3"4l . 
The second line of the redshift-distance relation (JBDJ) leads to 



(C30) 



„(2) 



£(i + e) njuj r dkk2p 

2ir{l + 2£) 2 A7 ?s Arjs, J a 



<E 

e=e~i,i+i 
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+ 1 



l- 



.(C31) 



Analogously the third line yields 

c (3) s £(i+J) i L_ 

* " 27r(l + 2£) 2 A?7s A775- 



dfc P CT (fc) 



€=£-i,£+r 



-(^-l) WA^-tJ (e(^_ 1)fc -»?s) e(^, fc -r?5') +0(r?l ifc -r? S ) e(^-_ M -^))]. (C32) 

The fourth line is composed of the two terms that contribute to 1Z given in Eq. (|55[) . Denoting them with superscripts 
(41), (42), and their correlation with (412) we obtain 



J4 i)„ 1(1+1) 1 1 



87r(l+2£) 2 Ar? s Ar? s 

,(42) s 1(1+1) 1 1_ 

8tt (1 + 2if Arjs Ar]s> 



-|dfc fc 2 p CT (fc)E [Pfa-vs) (% k -vs)ffi(% k )e(vzk-vs)e(vi, k -vs- 
1=1-1,1+1 



, (C33) 



dk P a (k) 



*E 



+ (»?!,* - (»?/_!,* " Vs>) - %?) (jll-x,k - , (C34) 



„(412) 



1 1 



dk k P a (k) 



8tt (1 + 2f) 2 A?7 S A77S' 

x E [(^,fc-w) (^-i.fe-^'J 1) h-iW* \ni-i,k) T k (m,k) (%k~vs) e (rn_ hk -vs' 
i=(-i,i+i 

- (vik-vs) (vik-vs') i(t + 1) ijTk (vik) T k (vik) (vik-ns) e (pik-m* 
+ (s <-> s") . 



(C35) 

Next, we compute the cross terms between the different lines of the distance-redshift relation ((59]) . We start with the 
second and third lines, 



„(23) 



t{l+l) Ug 1 1 



dk k P a (k) 



2tt (1 + 2lf A?7s A775. 

I h 7 k(vi.x,k) e ( r H-i,k- r is) &(m-i,k-vs' 



xE 

i=e-i.e+i 



k-ik T * (m-i,k) T k (vi,k) \ni-i,k - vs) (pit - Vs- 
^—j^k-iW* (m-i,k) T k (vi, h ) (%k - vs) {m-i,k - vs>) 



k T k [ T lLk - w ( w k - Vs 



+ {S^ S') . (C36) 



The second and fourth lines yield 
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x [t(ne,k -VS') 1 !- 



JiTk (jtt-i,k) T k \%k) 6 (vi-i,k - Vs) 6 (y l k - ns* 
is') ijTl (i7j, fc ) © (vi,k ~ Vs) (vik ~ Vs>) + 



- (* + 1) (m,k - vs>) ipl (v u ) e (vi, k - vs) e ( Vlk - V8 >)] + (s^ s') ., 

nil + 21? Atjs Ai 1s > J 

■)(* 



(C37) 



-Vs 
~ (% k ~ Vs- 

- (m-i,k - ? ?s 



Vs) (rii-x, k - Vs-) 



c (24 2) s ni+^j i dkkP (k) 

1 Air(l + 2£f At7s At7s' J * 

i=t-i,t+i 

l ) h-ih T * (m-i.k) Tk (vik) (^_i, fe - vs) e fe - tip) 

-= Tt: . hTi. ( Ht; . , ^ Tz. fr): , ^ B f?7j , — na\ B ^77; . , — 77c/ 



£ + 

') i —Y L h-\ I fi* (w-i,k) T k (ve >k ) © (yik - vs) © - vs-) 

(vik - VS') ^J-ljTk (v l>k ) f k ( %k ) B ( %k - m ) 6 ( %k - Vs ,) + (S o S>) , (C38) 



and, the third and fourth lines give 



,(341) 



t{l+l) 1 1 



- J dk k P a (k) 



„(342) 



4tt (1 + 2£) 2 A?7s At7 S / 

x J2 [(%k ~ Vs>) I (i- l) h-ih T * (vi-i,k) T k (vik) © (vi-i.k - Vs) (vi M - VS') 

e=i-i, m 

- (vik - vs') i(i+i) ipl ( Vik ) b ( %k - V s) (vik - vs' 

£(! + £) 1 1 



+ (S o S") , (C39) 



4tt (1 + 2£) 2 Arjs At? s 

2 



dfc P ff (fc) 



«-i,*fi 

- (^ 2 - l) (vik - VS') k-ihT k (m_ lik ) T k ( Vi k ) B (vi_ 1>k - vs) © (ve, k - Vs>) 



1 - 1 vi-ik - vs' h-ihn m k T k m k e ^ - vs © »z?_ lifc - ) 



^+ l) - »?/?') (r?i ife ) T fe (j7 i fe ) B (r} lk - rjs) B (r] lk - rjs') 



+ (S<+S'). (C40) 



To determine the correlation between the peculiar velocities and the shear on the constant time hypersurface 
we need to specify the gravitation theory. We choose GR by using the Einstein's equations (|5T|) . Correlating the term 
for the peculiar velocity of the source vg with the others, we find 



„( 12 ) 



1 



167rGa| (p s + p s ) 
+ (S o 5') 



i=i-i,t+i ls 



k(v)j'A kA v) 



1 Hs L \ Hp_ 2£(£+l) 



167rGa| (p s + Ps) V A? ?s / A?7 S / ^ (2/? + if 

x(^(, s + ^£)t 1 e(, s -w + ^ 



[¥ 4 P 






^A?7 S J 



T l/Ar ls ^s) 
£+1 

T i/AvAvs) ^-—=-6 fas - »7S') 



+ (5 o 5') , (C41) 



„(") 



•Hs 1 ^ 1 2^(^+1) 



167rGa|(/5 S +p s ) V A? 7s/ &VS'n {2£+lj 
+ (S o S") 



/■'fo 

dkk 5 P a (k)T k (vs) Ji( kA Vs) dn{vo-v)Tk{v)fi{kAn) 

7, ..... •'TJS' 



e=e-i,m 
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167rGa| 



T I/A VS I Vs- 



(Ps+Ps)\ Ar] S ) Ar] S , n (2£+l) 2 „ ^ 



A Vi 



,(141) 



X 



-i 
s 



(2£- 
-1) 6(t7s-?7s' 



e=e-i,e+i 
At] S 



Arf s 



Pa 



Ar 1S/ 



1 2£(£+l) 



8wGa 2 s (ps+Ps) V A Vs J A Vs' n (21 + if 

/rvo 
dk k 6 P a (k) T k ( Vs ) Yl H ( fcA ^) / d V (vo-v)(v-VS')T k fa) j'i (kAr,) 



e=£-i,e+i 



1- 



87rGa|(p s +P s ) V A? 7s / A??,?, ^ (2£ + l) 2 . 



1 2£(£ + l) 



1L V 



£—1,-6+1 L 



At?; 



-T; 



,(142) 



\ / 

l/Ar, s (.Vs)Ii (»?s-»7sO 9(r? s -r? s 0)] + (5 o 5') , 

1 2£(£+l) 



Hi 



87rGa| (ps+Ps) V Ar^/ A^ n (21+ 1 



/ dfc fc 4 P ff (fc) T fe ( Vs ) h ( fcA ^) / Mvo-v)(v-Vs')T k (r]) £ (fcA ?7 ) 

;„.„., J ils' 



e-i,e+i 



1 21(1 + 1) x ^ 
87rGa| (p s +P s ) V" A^y Ar^ ^ (2£+ 1) 2 ~_£^. 

+ (S <-> 5') . 









^A?7 S J 



+ (So 5') , (C42) 



+ (5 <-> 5') 



(C43) 



ji2 



+ (5 o 5') 
{ris-r]s>)<d(vs-Vs>', 



(C44) 



Appendix D: Limber approximation 770 — T) and At/ = 770 — 7] ) 

In this work we have used the Limber approxima- z^o j 
tion HH repeatedly. It approximates the integral of the / dr l T (v) k (kAr]) = —T (% jfe ) I t Q (r]^ k - r] S ) , (D3) 



no 



product of a spherical Bcssel function and a slowly vary- 
ing function (e.g., a power law) by f ° df] T ^ ^ / fcAr? ) 

J ri .q 



■1'2 



VS 

dxf(x)j e (x)^I e f(£)Q(x 2 -£)Q(£-x 1 ) (Dl) = \t (ru-i,h) h-i® (w-i,k - w) 

xi K 



r v T(n') 



- r T (Vt k) h^T-® (Vt k ~ VS) , (D4) 

for X2 > xi, where 9 denotes the Heaviside function de- k ' £ 

fined by p j , f (r)') 

, Vs (kA V ') 

( a; ) = { l>0 > ( D2 ) = ^f( mtk )I e e( m , k -r, s )@(v-Vi,k). (D5) 

and If = 1.58/ £ describes the area under the first peak / dr] / di) — — j t (kAr]') 

of the spherical Bcssel function je (x). This rather crude vs vs 



approximation considers the contribution under the first ~ \ rpi \t c\ /„ _ (r\a\ 

peak only, and it usually gives an overestimation, but K t 

never by more than a factor of 2 Of course, if the f no ^ , > , ? ) "H f ? ) T f? ) 3 l ^kAr]) 

function / varies heavily in the region of the first peak, J g ^ ^ ^' ^' (kAr]) 2 

£ — 1 < x < £ + 1 , the approximation cannot be used. ^ 

We are, in particular, interested in (note that Arj = - jj^Vt,k-Vs) IlH(Vt,k) T (W,k) 9 (r]e, k -Vs) , (D7) 
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mo £—1 

~k 

'Vs J vs \ 

= ^T(rit-i, k )It-iQ(ju-i,k-Vs) -^h(V£,k-Vs)T( Ve , k )Q( m , k - Vs ). (D10) 
t + 1 

fc2~ r IiQ ^t-k - Vs) , (D8) We haye uged the following p r0 p r i e ty f the spherical 

no Bessel functions [24j : 

ns 

- ^(VLk-Vs) ItTfak) e (VLk-Vs) , (D9) ft ( kA v) = jt-i (kArj) - —j £ (kArj) . (Dll) 

rno 

/ dr) (v-Vs)(vo-ri)T(v)fe(kAri) 
'vs 
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